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The total cross section for hadron production in e + e~ 
annihilation in the y/s region above the open-charm 
threshold was measured precisely by the Crystal Ball [l[ 
and BES [|J collaborations. However, the parameters of 
the J PC = 1 charmonium states obtained from fits 
to the inclusive cross section [!, 0] are poorly under- 
stood theoretically @. Since interference between dif- 
ferent resonant structures depends upon the specific fi- 
nal states, studies of exclusive cross sections for charmed 
meson pairs in this energy range are needed to clarify the 
situation. Recently, CLEO-c performed a scan over the 
energy range from 3.970 to 4.260 GeV and measured ex- 
clusive cross sections for DD, DD* and D*D* final states 
at twelve points with high accuracy Q . Belle has used a 
partial reconstruction technique to perform first measure- 
ments of the exclusive cross sections <r(e + e~ — > D ± Z)* =F ) 
and a(e+e- -> D*+D*~) for s/s near the D+D*~ and 
£)*+£)*- thresholds with initial-state radiation (ISR) 0. 
Recently Belle [|| has reported a measurement of the ex- 
clusive cross section for the process e + e~ — > D°D~ir + 
and the first observation of ^(4415) ->■ DD*.(2460) de- 
cay. 

In this paper we report measurements of the exclu- 
sive cross sections for the processes e + e~ — > D + D~ and 
e + e~ — > D°D° using ISR that are a continuation of our 
studies of the near-threshold exclusive open charm pro- 
duction. Recently several new charmonium-like states 
were observed in this mass range (F(4260) @, EH, 
y(4360), Y(4660) [HI, ^(4160) [h|) decaying to either 
open- or closed-charm final states. Our study provides 
further information on the dynamics of charm quarks 
at these center of mass energies. The data sample cor- 
responds to an integrated luminosity of 673 fb 1 col- 
lected with the Belle detector [l3| at the T(4S I ) resonance 
and nearby continuum at the KEKB asymmetric-energy 
e + er collider (l4l |. 

We select e + e~ — * DDjjsr signal events by recon- 
structing both the D and D mesons, where DD = D a D° 
or D + D~ . In general, the jisr is not required to be de- 



tected; its presence in the event is inferred from a peak at 
zero in the spectrum of the recoil mass against the DD 
system. The square of the recoil mass is defined as: 

M? ec (DD) = (E c . m . - E Dl5 f - p 2 Dls , (1) 

where -E c .m. is the initial e + e~ center-of-mass (cm.) en- 
ergy, E D jj and p D jj are the cm. energy and momen- 
tum of the DD combination, respectively. To suppress 
backgrounds we consider two cases: (1) the jisr is 
out of detector acceptance in which case the polar an- 
gle for the DD combination in the cm. frame is re- 
quired to be |cos(#£,-p)| > 0.9; (2) the fast jisr is 
within the detector acceptance (|cos(6' D ^)| < 0.9), in this 
case the jisr is required to be detected and the mass 
of the DDjisr combination is required to be greater 
than i? c . m . — 0.58 GcV/c 2 . To suppress background from 
e + e~ — > DD{n)'K r )isR processes we exclude events that 
contain additional charged tracks that are not used in 
the D or D reconstruction. 

We ensure that all charged tracks originate from the 
interaction point (IP) with the requirements dr < 2 cm 
and \dz\ < 4cm, where dr and \dz\ are the impact 
parameters perpendicular to and along the beam di- 
rection with respect to the IP. Charged kaons are re- 
quired to have a ratio of particle identification likeli- 
hoods, Vk = Ck/(£k + Ctt) [H], larger than 0.6. No 
identification requirements are applied for pion candi- 
dates. Kg candidates are reconstructed from 7r + 7r _ pairs 
with an invariant mass within 10MeV/c 2 of the nominal 
Kg mass. The distance between the two pion tracks at 
the K s vertex must be less than 1 cm, the transverse 
flight distance from the interaction point is required to 
be greater than 0.1cm, and the angle between the Kg 
momentum direction and the flight direction in the x — y 
plane should be smaller than 0.1 rad. The pion pair can- 
didates are refitted to the K s mass. Photons are recon- 
structed in the electromagnetic calorimeter as showers 
with energies greater than 50MeV that are not asso- 
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dated with charged tracks. Pairs of photons arc com- 
bined to form 7T° candidates. If the mass of a 77 pair lies 
within 15McV/c 2 of the nominal 7r° mass, the pair is fit 
with a 7T° mass constraint and considered as a ir° candi- 
date. D° candidates [l6| are reconstructed using five de- 



cay modes: K 7r + , K K + , K 



and 



K~n + n Q . D + candidates are reconstructed using the 
decay modes Kgir + and K~ir + Tr + . A ±15MeV/c 2 mass 
window is used for all modes except for K~ir~ir + ir + , 
where a ±10MeV/c 2 requirement is applied (~ 2.5 a in 
each case). To improve the momentum resolution of D 
meson candidates, final tracks are fitted to a common 
vertex and a D° or D + mass constraint is applied. The 
D candidates from a sideband region are refitted to the 
mass, corresponding to the center of the sideband region 
with the same width as the signal window. 

The distribution of M? CC (DD) after all the require- 
ments is shown in Fig. [T](a). A clear peak correspond- 
ing to the e + e~ — ► DD^/jsr process is evident around 
zero. The shoulder at positive values is due to e + e~ — » 
£)(*) D ( -*\n)Tr jisR events. To suppress the tail of such 
events we define a signal region by the requirement: 
\M? ec (DD)\ ^0.7(GeV/c 2 ) 2 . The polar angle distri- 
bution for DD after the requirements on M^ CC (DD) is 
shown in Fig.Q](c). The sharp peaking at ±1 is character- 
istic of ISR production and agrees with the Monte Carlo 
(MC) simulation. The mass distribution of the DDjjsr 
combinations where the jisr is detected is shown in 
Fig-ISO 3 ) after all the selection requirements. It peaks at 
the E c m . The asymmetric shape of the ^dTdjisr distri- 
bution is due to higher-order ISR processes. The M DO no 
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FIG. 2: The mass spectra of DD combinations (points with 
error bars): (a) D°D a ; (b) D + D~ . The total contribution 
from combinatorial background (1-2) is shown as the hatched 
histogram, the contribution from background (4) is shown as 
the (barely visible) solid histogram. 



The following sources of background are considered: 

(1) combinatorial D(D) mesons combined with a real 
D(D) coming from the signal or other processes; 



(2) both D and D are combinatorial; 
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FIG. 1: a) The observed distributions of (a) M? CC (DD); (b) 
M(DD'yisR) and (c)DD polar angles. Histograms show the 
normalized contributions from Md and Mjy sidebands. The 
selected signal windows are illustrated by vertical lines. 



and M D + D - spectra obtained after all the requirements 
are shown in Fig. [21 



(3) reflection from the processes e + e 

and e + e~ — ► DD*-fisR, followed by D 
with an extra 7r2,; c= in the final state; 



DDn^ iss -y IS R 



(4) reflection from the process e + e — > DD*^isr, fol- 



lowed by D*° 
final state; 



_D°7, with an extra soft 7 in the 



(5) a contribution from e + e — > DDir when an ener- 
getic 7T° is misidentified as a single Jisr- 

The contribution from background (1) is extracted us- 
ing Md and Mjy sidebands that arc four times as large 
as the signal region. These sidebands are shifted by 
30MeV/c 2 (20MeV/c 2 for the D° -> K-tt-tt + tt + mode) 
from the signal region to avoid signal over-subtraction. 
Background (2) is present in both the Mjj and M-^ side- 
bands and is, thus, subtracted twice. To account for this 
over-subtraction we use a 2-dimensional sideband region, 
where events are selected from both the Md and the Mjy 
sidebands. The total contribution of the combinatorial 
backgrounds (1-2) is shown in Fig. [5] as a hatched his- 
togram. Backgrounds (3-4) are suppressed by the tight 
requirement on M^ CQ {DD). The remaining contribu- 
tion from background (3) is estimated directly from the 
data by applying a similar full reconstruction method to 
the isospin-conjugate process e + e _ — > D° D~ ir^^jisR- 
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Here the requirement on absence of additional charged 
tracks in the event is relaxed. Since there is a charge 
imbalance in the D°D~ final state, only events with an 
extra missing 7r+ iss can contribute to the M 2 CC (DD) sig- 
nal window. To extract the level of background (3), the 
D°D~ mass spectrum is rescaled according to the ratio 
of D° and D~ reconstruction efficiencies and an isospin 
factor of 1/2. When this is done, the contribution from 
background (3) is found to be negligibly small. Uncer- 
tainties in this estimate are included in the systematic 
error. Background (4) contributes only to the D°D a 
final state. It is estimated using a MC simulation of 
e+e- -> D°D* j IS r, followedby D*° -► D°-f. To re- 
produce the shape of the D°D*° mass distribution we 
use the D D*+ cross section measured in our previous 
study Q . The contribution from background (4) is found 
to be small (shown in Fig.[2](a) as a solid histogram) and 
is subtracted from the D°D° mass spectrum. Uncertain- 
ties in this estimate are included in the systematic error. 
The contribution from background (5), determined from 



reconstructed e H 



DDir events in the data, is found 



to be negligibly small and taken into account in the sys- 
tematic error. 

The e + e~ — > DD cross sections are extracted from the 
L>°17 and D+D~ mass distributions [fl 



DD 



.TTN dN / dr 



i] tot dL/dm 



(2) 



where m = M D jj, dN/dm is the obtained mass spectra, 
while ?7tot is the total efficiency The factor dL /dm is the 
differential ISR luminosity 

dL/dm = — ( (2 - 2x + x 2 ) In - x 2 c) , (3) 
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FIG. 3: The exclusive cross sections for: (a) e + e — > D a D a ; 
(b)e + e _ -> D + D~\ (c) e+e~ -v ZXD. The dotted lines 
correspond to the ^(3770), ^(4040), ^(4160) and ^(4415) 
masses 12011. 



The systematic errors for the a (e+e — > Z?Z?) measure- 
ments are summarized in Table [I] The systematic errors 



TABLE I: Contributions to the systematic error in the cross 
sections, [%]. 



where x = 1 
luminosity and 



- m 2 /E 2 m , £ is the total integrated 
C = cos6>o, where 9q defines the po- 



lar angle range for the m the e + e~ cm. frame: 

9o < 1ISR < 180° — 8 . The total efficiency determined 
by MC simulation grows linearly with M D jj from 0.095% 
near threshold to 0.46% at 5GeV/c 2 for the D°D° and 
from 0.038% to 0_17% for the D+D- mode. The result- 
ing e+e~ -> D°D°, e+e" -> D+D~ and e+e" -> L>D 
exclusive cross sections, averaged over the bin width, 
are shown in Fig. [3] with statistical uncertainties only. 
Since the bin width in the cross section distributions is 
much larger than the resolution (which is ~ 3MeV/c 2 
at threshold and ~ 5McV/c 2 at ~ 5GeV/c 2 ), no 

correction for resolution is applied. 

We calculate the cross section ratio 
a(e+e- -> D+D- )ja(e+e~ -> D°W) for the M m bin 
(3.76 — 3.78) GeV/c 2 corresponding to M D -^ m M^,(3 77 o) 
to be (0.72 ± 0.16 ± 0.06). This value is in agreement 
within errors with CLEO-c and BES measure- 
ments. The ratio a(e+e~ -> D+ D~) / o{e+ e~ -> D°D°) 
integrated over the M D jj range from 3.8 to 5.0GeV/c 2 
is found to be (1.15 ± 0.13 ± 0.10) and is consistent with 
unity. 



Source 




D+D- 


DD 


Background subtraction 


±4 


±3 


±3 


Reconstruction 


±7 


±6 


±7 


Cross section calculation 


±5 


±5 


±5 


B(D) 


±4 


±6 


±5 


Kaon identification 


±2 


±2 


±2 


Total 


±10 


±10 


±10 



associated with the background (1-2) subtraction are es- 
timated to be 2% due to the uncertainty in the scaling 
factors for the sideband subtractions. This systematic 
error is estimated using fits to the Mq and M-p distri- 
butions with different signal and background parameter- 
izations. Uncertainties in backgrounds (3-5) are conser- 
vatively estimated to be smaller than 2% of the signal in 
the case of D°D°; these two sources are added linearly to 
give 4% in total. In the D+D~ case, backgrounds (3-5) 
are estimated using the data and only the uncertainty 
in the scaling factor for the subtracted distributions is 
taken into account. A second source of systematic error 
comes from the uncertainties in track and photon recon- 
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struction efficiencies, which are 1% per track, 1.5% per 
photon and 5% per Kg, respectively. The systematic er- 
ror ascribed to the cross section calculation is estimated 
to be 5% and includes the error on the differential ISR 
luminosity and the error from the efficiency fit. Other 
contributions come from the uncertainty in the identifi- 
cation efficiency and the absolute D° and D + branching 
fractions [2(| ■ The total systematic uncertainties are 10% 
and comparable to the statistical errors in the differen- 
tial cross section around the -0(3770) peak; for the other 
M D j=j ranges statistical errors dominate. 

In summary, we report measurements of e + e~ — > 
D°D° and e + e~ — > D + D~ exclusive cross sections for 
J s near the D°D° and D+D~ thresholds with initial- 
state radiation. The observed e + e~ — > DD exclusive 
cross sections are consistent with recent BaBar measure- 
ments [2~l| and arc in qualitative agreement with the 
coupled-channel model predictions of Ref. [23]. This in- 
cludes the peak at 3.9GeV/c 2 that is seen both in Belle 
and BaBar cross section spectra. 
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